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ABSTRACT 



A generic prediction of the single-degenerate model for Type la supemovae (SNe la) is 
that a significant amount of material will be stripped from the donor star (~0.5 Mq for a giant 
donor and ~0. 15 M0 for a main sequence donor) by the supernova ejecta. This material, excited 
by gamma-rays from radioactive decay, would then produce relatively narrow (<1000kms-i) 
emission features observable once the supernova enters the nebular phase. Such emission has 
never been detected, which already provides strong constraints on Type la progenitor models. In 
this Letter we report the deepest limit yet on the presence of Ha emission originating from the 
stripped hydrogen in the nebular spectrum of a Type la supernova obtained using a high signal- 
to-noise spectrum of the nearby normal SN la 201 Ife 274 days after B-band maximum light 
with the Large Binocular Telescope's Multi-Object Double Spectrograph. We put a conservative 
upper limit on the Ha flux of 3.14 x 10~^^ er g s~^ cm~^, vyhich c orresponds to a lumi nosity of 
1.57 X 10^^ erg s~^ Assuming the models of iMattila et al.l (|2005h and the methods of Leonard 
(120070 . this translates into an upper limit of < 0.001 Mp , of str ipped material, which is an order 
of magnitude stronger than previous limits by [Leonard (l2007h . SN 201 Ife was a typical Type 
la supernova, special only in its proximity, and we argue that lack of hydrogen emission in its 
nebular spectrum adds yet another strong constraint on the single degenerate class of models 
for SNe la. 
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1. Introduction 



Despite the fact that SNe la were used to discover the accelerating universe (iRiess et al.lll998l : lPerlmutter et al. 



19991) ■ the physical nature of th eir progenitor systems remains theoretically ambiguous and observationally 
elusive (for a review see Wang & Hanll2012l ). It is commonly accepted that SNe la result from the ther- 
monuclear explosion of a carbon-oxygen white dwarf (WD) in a close binary system, but the nature of 
the binary companion and the sequence of events leading to the SN explosion are still uncertain. There 
are two dominant models: the double degenerate (DD) scenario, in which the companion is also a WD 
dlben & Tutukovl 1 19841 : IWebbinld 1 19841) . and the single degenerate (SD) scenario, in which the compan- 
ion i s a non-degenerate ob j ect such as a m ain sequence (MS) star, a red giant (RG), a sub-giant or a He 



star (IWhelan&Ibenlll973l : 



Nomoto 



1982h . There are additional minor variations on these basic models 



of the p rogenitor system for both the SD channel (e.g . ljusthami l201ll : IWheelen l2012h and the DD chan- 
nel (e.g. lThompsonll201 ll : IShappee & Thompsonll2012h . These uncertainties about the progenitor systems 
for SNe la rem ain a substantial problern for understanding the systematic errors in using SNe la to study 
cosmology (e.g. IWood-Vasey et allbOOTT) . 



One of the observational signatur es of the SD model is that the material from the companion should be 
stripped when struck by the SN ejecta (|Wheeler et al.lll975h. leading to both imrnediate signatures from the 
impact of the SN ejecta on the co mpanion (e.g. Marietta et al.ll200Cll : iMeng et al.ll2007l : IPakrnor et al.ll2008l: 



Pan et al 



and 



20121 



3; 



Liu et al 



Liu et al 



2012!) an d majo r changes in the companion's future evoluti on Jpodsiadlowskil 



2003l:lshappee et al.)l2012l : |Pan et al]|2012al) . Recently, the hydrodynamic simulations of |Pan et all (l2012bl) 



(120 121) have shown that ~ 0.1-0.2 Mq of solar-metallicity material is expected to be re- 
moved from MS companions by the impact of the SN ejecta. These hydrodynamic simulations show that 
this material will be embedded in low-velocity supernova debris with a characteristic velocity of <1000 
km s~' . However, the line profile s from the stripp ed material are somewhat uncertain because this material 
will be asymmetrically stripped (|Liu et al.ll2012h . and the orientation of the binary relative to our line of 
sight, at the time of explosion, is unknown. This material will be hidden in early-time spectra by higher 
veloc ity, optically thick, iron-rich ejecta, but will then appear in late-time, nebular phase spectra (> 250 
days: lMattila et al.ll2005h as the higher velocity ejecta become optically thin. 



Only a handful of nebular phase, high signal-to-noise (S/N) SNe la s pectra have been pu blish ed in the 
literatu r e, with the strongest limits on late-time hydrogen flux coming from lMattila et all (l2005h and lLeonard 
(|2007h . iMattila et all (|2005h obtained late-time low-resolution spectra of SN 2001el, modeled the emission 
from solar-metallicity material stripped from a non-de generate compan ion, and used both to place an upper 
limit of <0.03 Mq on the presence of this material. Leonard! (120070 obtained deep, medium-resolution, 
nebular-phase, multi-epoch spectra of SN 2005am and SN 2005cf at distances ~37 Mpc and ~32 Mpc, 
respectively. Using the approach of iMattila et all (l2005h . he placed a <0.01 Mq upper limit on the presence 
of any low velocity solar-abundance material in both SN. These observations seem to firmly rule out main 
sequence companions for these three SNe. 



However, a more exotic scenario which might evade these constraints has since been proposed by 
JusthamI (|2011r) . In this model a WD is spun up by the mass it accretes from its non-degenerate companion. 
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allowing the WD to remain stable above the Chandrasekhar mass and giving the companion time to evolve 
and contract before the SN explosion. This leaves a smaller and more tightly bound companion at the time 
of explosion, reducing the amount of stripped material. It therefore is of significant interest to put even 
stronger constraints on the amount of stripped material for other SNe la. 



At a mere 6.4 Mpc (|Shappee & Stanekll201 ih the "plain vanilla" Type la SN (as described by [Wheeler 



2OI2I) is an ideal target for obtaining improved c onstraints. SN 2 01 Ife was discovered less than one day 
after explosion by the Palomar Transient Factory (ILaw et al.ll2009l) in the Pinwheel Galaxy (Ml 01), a well- 
studied, nearby face-on spiral. This SN is the nearest la in the last 25 years and the early detection and 
announcement allowed extensive multi-wavelength follow-up observations. Additionally, SN 2 life is only 



slightl y red dened and is surro unded by a relatively "clean" environment (IPatat et al.l 120111) . iBrown et al. 



(120 12I ) and Ib loom et all (120121) used the very early-time UV and optical observations to constrain the shock 
heating that would occur as the SN ejecta colhde with the companion, ruling out giants and main sequence 
secondaries more rnassive than the Sun. These constraints appear to rule out SD models; however, these 
studies both rely on lKasenlJ2010 l). which assumes a Roche-lobe overflowing secondary with a ty pical stellar 



structu re and depends on the orientation of the binary relative to our line of sight. However, iMeng et al. 



(|2007h pointed out that the secondary will be left smaller and more compact than a typical main sequence 
star after binary evolution, which would produce a more subtle photometric signal at even earlier times 
which could have been missed by early-time observations. 

In this study we place the strongest limits yet on the presence of Ha emission in the nebular spec- 
trum of a Type la supernova (SN la). We obtained a very high S/N spectrum of SN la 201 Ife 274 days after 
maximum B-band light using the Multi-Object Double Spectro graph on the 8.4 -meter Large Binocular Tele- 
scope. We generally followed the analysis procedures used by Leonard! (|2007n . so our constraints should be 
directly comparable to earlier results. In Section |2] we describe our observations, in Section [3] we derive our 
Ha limits, and we summarize our findings in Section IH 



2. Observations 

We obtained five high-S/N spectra of SN 201 Ife from 73 to 274 days after maximum B-band light 



5ing the first of the Multi-Object Double Spectrographs (MODS;|Poggeetal ]l2O10h on the Large Binocular 



usm 

Telescope (LBT). The two channels of MODS allowed us to obtain wide spectral coverage (3200- lOOOOA) 
in a single exposure, which will make these observations useful for many future studies. Details of the 
observations and the flux standards used are listed in Table [T] All observations were taken through a 1" wide 
slit. To perform basic CCD reductions on our spectra we followed the "MODS Basic CCD Reduction with 



modsCCDRed" manual^. We then performed cosmic ray rejection using L.A.Cosmic ( van Dokkum 2001 ) 



and combined each channel's 2D spectra for each epoch. Next, we extracted the ID sky-subtracted spectra 



^modsCCDRed manual is available here |http: / /www, astronomy, ohio- state . edu/MODS /Manual s/MODSCCDRed.p'df| 
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using the apall task in Each spectrum was then waveleng th and flux caUb rated. We did not attempt 

to correct small-scale telluric absorption bands as performed by [Leonard (120071) . This will not adversely 
affect our results because the telluric absorptions occurs on wavelength scales smaller than the Ha emission 
expected from the SN. This also avoids adding any noise from telluric spectra to our SN spectra. 

Because our spectra were taken under non-photometric conditions with a relatively narrow (1") slit, it is 
necessary to scale the fluxes of our spe ctra to place them on an absolute flux scale. To do this we performed 
aperture photometry on our Sloan r' (IFukugita et al.l Il996h acquisition images using the IRAF package 
apphot. Because the SN was saturated in the acquisition image of the first epoch we must treat that spectrum 
separately. For the remaining epochs we used the four b rightest stars in t he image to put our pho tometry on 
the Sloan Digital Sky Survey (SPSS: I York et allboOOl) Data Release 7 dAbazajian et allboogf) magnitude 
system. Our r'-band magnitudes are reported in Table [T] We then scale the spectrum so that its synthetic 
r'-band photometry matches the computed /-ban d aperture photometry. F or the first epoch we scaled the 
spectrum such that a synthetic /?-band magnitude (|Bessell & Murphvll2012h matches the j?-band magn itude 
for its epoch found by interpolating the combined /?-band light curves of [Richmond & SmithI (|2012l ') and 
Munari et al. ( 2012 ). To convert from the AB magnitude system (OkQ 19741) of our synthetic photometry 



to the Vega magnitude s ystem reported inlRichmond & SmithI (|2012r) and lMunari et al.l (|2012r) . we use the 
conversion presented in Blanton & RoweisI (2007). The factor by which each spectrum was multiplied is 
reported in Table [T] From the last four epochs' spectra, whi ch are calibrated from the acqu isition image 



but are also covered by the combined /?-band light curves of Richmond & SmithI (120 12h and iMunari et al. 



(l2012h . we estimate that our absolute flux calibration in the j?-band (where Ha i s loca ted) is accurate to 



10% o r better. In Figure [T] we show the BVRI light curves of [Richmond & SmithI (l2012r) and iMunari et al. 



(|2012|) as compared to our synthetic BVRI photometry from the spectra for illustrative purposes. Figure 
[2] presents our calibrated spectra with ob vious telluric fe atures marked. Finally, we compare our nebular 
phase spectrum to the spectra presented in lLeonardI (|2007|), which provided the previous best limits on Ha. 



'iRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with the National Science Foundation. 



Table 1. Observations 







HJD 


P.A. 


Par. P.A. 






Flu.x 


Seeing 


Exposure 




Scale 


UT Date 


Day 


-2,400,000 


(dcg) 


(deg) 


Airmass 


standards 


(aresec) 


(s) 


(mag) 


factor 


2011 Nov 23.52 


73.02 


55889.02 


-78.0 


-84.3--85.7 


1.60- 


1.64 


FeigellO,G191-B2B 


0.8/0.9 


360.0/360.0 




1.07 


2012 Jan 2.54 


113.05 


55929.05 


-124.4 


-134--139 


1.11 - 


1.12 


Feige34 


0.8/0.9 


900.0/900.0 


14.26 ±0.03 


1.93 


2012 Mar 24.47 


194.97 


56010.97 


-68.5 


110-124 


1.15 - 


1.24 


Feige67 


0.5 / 0.7 


2880.0/2160.0 


16.35 ±0.01 


1.77 


2012 Apr 27.27 


228.77 


56044.77 


-2.0 


-159-177 


1.07- 


1.08 


HZ44, BD+33d2642 


0.9/1.1 


3600.0/3600.1 


17.06 ±0.01 


2.20 


2012 Jun 12.16 


274.66 


56090.66 


-32.7 


-177-129 


1.07- 


1.13 


HZ44 


0.7/0.7 


7200.0 / 7200.0 


17.95 ±0.01 


1.26 



No te. — Observational and derived properties of our time series spectra of SN 201 Ife. Days since maximum B brightness assume JD rgrnax = 2455816.0 zb 0.3 [Richmond & Smitl] 
I2OI2I) . P.A. is the position angle of the spectrograph slit. Par. P.A. and airmass give the range of parallactic angles and airmasses at the start of each separate observation, respectively. 
Standai^d stars were observed on the same night as science observations and, when multiple standai^ds were available, the computed response functions from each standai^d were averaged. 
Seeing gives the FWHM for the red / blue channels of the spatial profile in the combined 2d spectra. Exposure times are given for the red / blue channels sepai^ately. r magnitude is derived 
from the acquisition images. Each flux calibrated spectra was multiplied by a scale factor derived in the R- or r -band to place it on an absolute flux scale. 

'^SN 201 Ife was saturated in this epoch's acquisition images, so no reliable photometry was obtained. 
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As can be seen in Figure [2l our spectrum has a substantially higher S/N than the previous studies, 
observations allow us to place even stronger limit on Ha emission as we discuss in ^ 
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Fig. 1. — SN 201 Ife hght curv es relative to the time of maximum B-band brightness (fsmax = 2455816.0 ± 
0.3. Richmond & SmithI l2012h. Open diarnonds and open squares show the BVRI photometry from 
Richmond & Smithr ( 2012 ) and Munari et al. (2012), respectively. Our synthetic BVRI photometry are rep- 
resent ed by filled stars. The first epoch is calibrated to an absolute flux scale using the ^ Richmond & Smith 
(120 12h /?-band light curve so good agreement is expected. The rest of our spectra are calibrated using our r'- 
band acquisition images and SDSS DRV photometry. From the second epoch we estimate that our absolute 
flux calibration is accurate to better then 10% in the R-hand. 



3. The Search for Hydrogen 

In our search for Ha emission in the spectrum of SN 201 Ife, we closely follow the methods presented 



m lLeonardl (120071) . We define a continuum by smoothing our spectra on scales large compared to expected 
Ha feature, then subtract off this continuum and examine the residuals. We searc hed for Ha emission withi n 

1991'). 



±1000 kms"^ (±22 A) about Ha at theredshift of MlOl, 0.000804±0.000007 (ide Vaucouleurs et al. 



We smoothed the spectrum with a second-order Savitsky-Golay smoothing polynomial (IPress et al. 



1992) 

with a width of 60 A. Our data requires this smaller smoothing scale than that employed by [Leonard! (12007.) 
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4000 6000 8000 10000 
Rest Wavelength (A) 

Fig. 2. — Rest-frame absolute-flux-calibrated spectra of SN 201 Ife. Vertical dashed gray lines mark obvious 
telluric features. Spectra are labeled by days since maximum B brightness. 
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Fig. 3. — Comparison b etween our nebu lar phase spectrum of SN 201 Ife and the spectra of SN 2005cf and 
SN 2005am presented in [Leonard (|2007l ) which previously provide d the best l imits o n Ha emission from the 
stripped material. Our spectra are binned to match the resolution of lLeonardl (|2007h . The rest wavelength of 
Ha is indicated by the vertical red line and the shaded gray region shows where hydrogen emission would 
be expected ( ±1000 kms~' = ±22 A about Ha). Our spectrum has a substantially higher S/N due to the 
proximity of SN 201 Ife, allowing us to place a substantially stronger limit on Ha emission in Section [3l 
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for a good continuum fit. It is, however, still significantly larger than the expected velocity width of any Ha 
feature. Additionally, we smoothed over the telluric feature at 6510-6525 A before smoothing the rest of 
the spectrum, because the feature was affecting the continuum determination near Ha. 

Our nebular phase spectrum and continuum fit are shown in the top panel of Figured in the vicinity of 
Ha and are binned to their approximate spectral resolution. 

We then subtracted the continuum from the binned spectrum and examined the residuals, shown in 
the bottom panel of Figure |4l for narrow Ha emission. There is no evidence for any Ha emission in the 
spectrum. The closest possible emission is the hump in the smoothed continuum at '--'6575-6595 A which 
is discussed in Section [l!2l More broadly, we found no narrow emission lines at any wavelength, including 
regions around H/3, [01] AA6300,6364, [Oil] AA73I9,7330, [O III] AA4959,5007, [Call] AA729I,7324 , 
and [Nil] AA6548,6583. 



3.1. Statistical Limit 



Following [Leonard & Filippenkd (120011) and [Leonard (120071) . we compute a 3a upper bound on the 
equivalent width as 



Wx(3a) = 3AA A/ThWATTITb 
= 3A/ywWAX, 



(1) 



where A A is the width of a resolution element (in A), AI is the la root-mean-square fluctuation of the flux 
around a normalized continuum level, Wiine is the full-width at half-maximum (FWHM) of the expected 
spectral feature (in A), B is the number of bins per resolution element in the spectrum, and AX is the bin 
size of the spectrum (in A). For the spectrum shown in in Figure IH AI = 0.0024 A and AX = 4 Aleading to 
Wx{3a) = 0.067 A for W^e = 22 A. 

To translate this equivalent width constraint int o a constraint on the a mount of material st ripped from a 
non-degenerate companion we follow the analysis of lMattila et all (|2005l ). iMattila et all (l2005h estimate that 
0.5 M0 of solar-abundance material 380 days after explosion produces an Ha feature with peak luminosity 
of '--'3.36 X 10^^ ergs s ~' A~^ . Accounting f or the distance to MlOl and Galactic extinction towards MIOl of 
E{B-V) = 0.009 mag (|Schlegel etal.lll998h and assummg /?y = 3.1, the expected Ha peak flux from 0.5 M© 



of stripped material in SN 201 Ife is 6.71 x 10"'^ erg s"' cm~^ A~^ If we assume the line is Gaussian with 
FWHM 22 A, then this feature would have an equivalent width of Wa(0.05 Mq) = 5.8 A. Scaling linearly 
from the equivalent width of the Ha emission line to the amount of stripped material, we place an upper 
limit on the amount of solar-abundance material in SN 201 Ife of 5.8 x 10"^ Mq. 
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Fig. 4. — Nebular phase spectrum of SN 201 Ife illustrating our conservative limit of 3.14 x 10~'^ erg 
s~^ cm~^ on the flux of Ha emission. The rest wavelength of Ha is indicated by the vertical red line 
and the shaded gray region shows where hydrogen emis sion would be expected ( ±1000 kms~^ = ±22 A 
about Ha). Adopting the models of iMattila et all (l2005h . these limits translate into a <0.001 Umit on 
the amount of solar-abundance material stripped from the companion. Vertical dashed gray line marks a 
large telluric feature. Top Panel: SN spectrum binned to the approximate spectral resolution (4.0 A; black 
solid); smoothed continuum (solid red); smoothed continuum with Ha limit added (dashed red); smoothed 
continuum with Ha limit subtracted to show what the underlying continuum would have to be to get the 
observed spectrum (dashed blue); and smoothed continuum with Ha limit subtracted, assuming the velocity 
distribution discussed in Section 13.21 (dotted blue). Middle Panel: The location of telluric water vapor 
absorption lines illustrated by the ESO SM-01 Sky Model Mode Version 1.3.1. Bottom Panel: SN spectrum 
with smoothed continuum subtracted (solid black) as compared to the Ha limit (dashed red). Horizontal 
solid red line marks zero. 
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3.2. Conservative Limit 



Unlike the iLeonardI (120071) study, the main uncertainty in our Ha limit arises from the continuum 
determination and the Ha line profile rather then the photon noise of our spectrum. For example there is 
a small amplitude feature in the observed spectra at 6575 - 6595 A. To test if this spectral feature could be 
attributed to material asy mmetrically stripp ed from a companion we took the velocity distribution of stripped 
mater ial from Figure 9 of lLiu et alj (120121 ) and assumed it was narrowly distributed along a line (see Figure 
11 of iLiu et alj|2012l) . We found that a total flux of 2.20 x 10"^' erg s"^ cm~^ and an angle between the 



stripped material and our line of sight of 30.0° best reproduces the feature in our smoothed continuum. To 
be conservative, we take a less stringent limit on the total Ha flux of 3.14 x 10"^' erg s"^ cm"^, which is 
weaker than our statistical limit and too laige to explain the spectial feature between 6575-6595 A shown in 
FigureHl This limit corresponds to an equivalent width of 0.12 A, a to tal Ha luminosity of 1.57 x 10^^ erg 



s ^ and 0.001 Mq of solar-abundance material adopting the models of iMattila et al.l (120051) . We emphasize 



that whether this feature is Ha emission or dumpiness in the underlaying SN continuum, our conservative 
flux limit holds. 



4. Summary 



We obtained five deep, medium resolution, spectra of the nearby "plain vanilla" SN la 201 Ife with 
LBT/MODS from 73 - 274 days after maximum B-band light. With the last nebular phase spectrum we 
place the deepest flux limits on narrow Ha emission yet for a SN la. Determining the late-time Ha emission 



inSNe la spectra requires difficult radiative transfer calculations, but adopting the models of iMattila et al. 



(120051 ). our limit translates into an upper limit on th e mass of sola r-abundance material of <0.001 Mq, 
an order of a r nagnitude smaller th en previous limits (ILeonardI 120071) . However, two important theoretical 
questions from IMattila et all (120051) need to be confirmed by additional modeling. First, the opacity of the 
high-velocity iron-rich ejecta should be recomputed to confirm when the SN ejecta becomes optically thin. 
Second, the expected excitation for Ha e mission by gamma-r ays should be modeled in detail. These issues 
notwithstanding, our limit based on the IMattila et al.l (120050 models poses a significant challenge to SD 
models, including more exotic variations proposed to evade the previous constraints. 
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